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Abstract 

We examine a potential role of the neutrino deuteron reactions in the mechanism of supernova 
explosion by evaluating the energy transfer cross section for the neutrino heating. We calculate 
the energy loss rate due to the neutrino absorptions through the charged-current process as well as 
the neutrino scattering through the neutral-current process. In so doing, we adopt a detailed eval- 
uation of cross sections for the neutrino deuteron reactions with the phenomenological Lagrangian 
approach. We find the energy transfer cross section for the deuteron is larger than those for ^H, 
^He and ^He for neutrino temperatures {Tu ~ 4 MeV) relevant to supernova core. Because of the 
low energy threshold for the deuteron breakup, the energy transfer rate rapidly increases from low 
temperature, 1 MeV. This suggests that the neutrino deuteron reactions may contribute ef- 

fectively to the heating mechanism during the dissociation of irons into light elements and nucleons 
in the shocked material of supernova core. 

PACS numbers: 21.45.Bc, 23.40.Bw, 25.30.Pt, 26.50. +x, 97.60.Bw 
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I. INTRODUCTION 



The reactions between neutrinos and nuclei play important roles in various phenom- 
ena in physics and astrophysics. The neutrino-nucleus interactions are crucial in physics 
of core-collapse supernovae through neutrino scattering, absorption and emission. Reliable 
evaluation of the rates of neutrino-nucleus interactions is necessary for the terrestrial neu- 
trino detectors such as Super-Kamiokande and SNO [l|, 0, H, 0, El, 0], for predicting the 
amount of elements produced in the nucleosynthesis through neutrino-processes 0, HI, 0], 
and for understanding the mechanism of supernova explosion [lol. 11] after the core bounce 
of the gravitational collapse of massive stars. 

In the neutrino heating mechanism, which is one of the key issues for the successful 
supernova explosion, the neutrino-nucleus interactions may contribute to the energy depo- 



sition in addition to the neutrino-nucleon interactions |lll . |12| . In most of recent studies 
on supernovae, the shock wave launched by the core bounce stalls on the way and needs 
extra assistance to revive the outward propagation 13, 3]- Although the absorption of 



neutrino on the nucleon is the main mechanism of the energy deposition, nucleus can be an 
extra agent of the energy deposition through the absorption or scattering. An advantage of 
nuclei as targets is that neutrinos interact through the neutral-current (NC) as well as the 
charged- current (CC). Since the neutrinos of /i- and r-types have relatively higher average 
energies than those of electron- type, they can have larger cross sections than those from the 
CC reactions for the low-energy electron-type neutrinos. 

The neutrino heating through the neutrino-nucleus interactions has been proposed by 
Haxton and the average energy transfer cross section has been evaluated for the representa- 



tive case of '^He and ^^Fe ll|. Since the matter from outer layer free-falls down to the shock 
wave and is dissociated into ^He and nucleons due to the shock heating, those two represen- 
tative species are potential targets of neutrinos streaming from the central core. Although 
its influence on the supernova dynamics has been studied by numerical simulations under 
the spherical symmetry with those reaction rates l2], it was reported that their effect is 



too small to affect significantly the dynamics. In a recent study on the standing accretion 
shock- wave instability (SASI), Ohnishi et al. [l^ have explored how the inelastic scattering 
of neutrinos with ^He influence the outcome of shock dynamics by performing numerical 
simulations in two dimensions. An effect on the growth of instability toward the explosion 
is found only for the critical condition in which the dynamics of shock wave is marginal 
to the explosion and exploratory cases with enhanced reaction rates. The effect depends 
crucially on the cross section of the neutrino-^He scattering and detailed evaluation based 
on theoretical models have been done recently j^, [l6| . 

Recently, several authors pointed out the appearance of light elements (A < 3) other 
than ^He in supernova environment, and its potential role on the supernovae through the 
neutrino interactions 17. Hsl. Il9|. O'Connor et al. [l3l studied the abundance of tritons and 
^He in dense matter and evaluated the neutrino reaction rates. The A=3 nuclei can be more 
abundant than ^He and nucleons, and its energy transfer rate can be larger than that of ^He 



for certain conditions of proton fraction, density and temperature. Arcones et al. [19| used 
a statistical model and studied the abundance of light elements in the neutrino-driven flows 
after the explosion. They demonstrated that the neutrino reactions on the light elements 
can have an influence on the spectra of neutrinos emitted from the surface of proto-neutron 
stars. 

So far, attention has been paid to the neutrino reactions on A=3, 4 nuclei only. Sumiyoshi 
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and Ropke [181] studied the composition of light elements in supernovae by including the 
deuteron, A=3, 4 nuclei and others in many-body calculations (See also [l9|), and found 
that the deuteron can also appear abundantly in some regions of supernova environment. 
Since the cross section of the neutrino-deuteron reaction is much larger than those for ^H, 
^He and ^He, the energy deposition through the deuteron may significantly contribute to 
the neutrino heating during the dissociation of iron nuclei into light elements and nucleons. 

In this paper, we examine the possible role of the deuteron by evaluating the neutrino 
heating rates through the neutrino-deuteron (ud) interactions. For the estimation of the 
ud reaction cross sections, we adopt the standard nuclear physics approach (SNPA)[3l, 3]. 
In SNPA, we consider the nuclear current consisting of one-body impulse terms and two- 
body exchange-current terms, and evaluate their matrix element with the nuclear wave 
functions generated with a high-precision NN potential. The reliability of our SNPA has 
been extensively tested by comparisons of calculated observables with available data of 
photo-processes, and also from a good agreement with the results of the effective field theory 
approach 0, 0] for the solar neutrino energy. We calculate the average energy transfer rates 
(and the heating rate) by using the Fermi distribution of the neutrino energy, and examine 
an effect of the ud reaction on the supernova dynamics in the heating region. We have in 
our mind here the typical situation realized in the region between the stalled shock wave 
and the proto-neutron star surface, where neutrinos are emitted with spectra close to the 
one under the equilibrium parametrized by the local temperature. We consider the neutrino 
reactions on the deuteron through both the CC and NC, covering the deuteron breakup and 
scattering. 

We show that the energy transfer rates through the ud reactions are much larger than 
those of ^H, ^He and ^He for the typical range of the neutrino temperature, especially at 
low temperature (energy) because of the lower breakup threshold. This suggests that the 
deuteron can contribute, as an extra agent, to the neutrino heating in supernova dynamics. 
The current finding urges ones to examine its effects on the supernova explosion by imple- 
menting the ud reactions as well as a detailed treatment of the mixture of light elements in 
supernova simulations. 



II. FORMULATION 

We are concerned with the charged- current (CC) and neutral-current (NC) v[i')d reac- 
tions listed in the following equations. 



Ve + d^e~+p + p [uCC] (1) 

Ue + d-^e^ + n + n [uCC] (2) 

u + d^u + p + n [uNC] (3) 

u + d^u+p + n [uNC] (4) 

u/u + d^u/u + d [u/uscatt] (5) 



In the supernova environment, the CC reactions act as an absorber of the neutrino, and the 
neutrino energy is deposited to the rest of compositions of supernova matter. For the NC 
reactions, the energy difference between the incoming and outgoing neutrinos is the energy 
transfer to the matter. 

The interaction Hamiltonian for semileptonic weak processes is given by the product of 
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the hadron current (Ja) and the lepton current (L^) as 



^cc ^ I dx[J^^{x)L^^'\x) + h. c] , (6) 



for the CC process and 



NC 



^ldx[jrix)L^^^\x)+h.c.], (7) 



for the NC process. Here G'p is the Fermi couphng constant, and Vud is CKM matrix element. 
For the weak couphng constant we adopt G'p = 1.1803 x 10~^GeV"^ 0|. The CKM matrix 
element is taken to be Vud = 0.9740. 
The hadronic CC is written as 

J^^ix) = V^^^\x) + Al^^\x), (8) 

where V\ and Ax denote the vector and axial-vector currents, respectively. The J^^*^ for 
the z//i^-reaction denotes ± iJ"^, where J* is the ith component of the isovector current. 
The hadronic NC is given by 

J^^{x) = (1 -2sm^ew)V^{x) + Al{x) - 2 sin^ OwV^ (x) , (9) 

where 6w is the Weinberg angle and is the isoscalar part of the vector current. The 
lepton currents, L^^'^ and L^*-^'^, are well known. 

The nuclear current consists of one-nucleon impulse approximation (lA) terms and two- 
body meson exchange-current (EXC) terms. The lA current is determined by the single- 
nucleon matrix element of Ja with the standard parametrization. For the axial vector 
EXC {Aexc)i we consider the pion-pair current, rho-pair current, pion- and rho-exchange 



A currents, n — p current following Refs. (20|, l21[. The strength of Apxc is adjusted 



to reproduce the experimental value of the triton beta decay rate. Regarding the vector 
EXC, we take into account the pion-pair, pionic and pion- and rho-exchange A currents. 
As discussed in jsf, the model of the nuclear vector current leads to np — > o?7 total cross 
sections that agree very well with the experimental values. The expressions of lA, EXC and 
the coupling constant used in this work are given in Ref . [1, 0] . 

The cross sections for the u/u^k) + d{P) l{k') + Ni{p[) + -/V2(P2) in the laboratory 
system are calculated following the standard procedure. We obtain the cross section for the 
CC reaction as 

da = E ^'''^''^l^^^/'^'^^'f'' E',) \lY.'\^dk'dp[dp',, (10) 

where we have included the Fermi function F{Z, E[) to take into account the Coulomb 
interaction between the electron and the nucleons in the final state. The cross section for 
the NC reaction is written as 

= E ^'''^^ ""(2 ^5^^ " ""^ ^rfx^l-dk'dpW,. (11) 
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We have used the matrix elements and defined as 

jA = <NN{P')\.UOMP)>, (12) 
h = <lik')\Lx\iy/uik)> . (13) 

The cross section for the neutrino-deuteron elastic scattering {^/^{k) + d{P = 0) — >■ 
v/vik') + d{P')) is written as 

^ = '^M£G'^4smHw[Acos'eL/2 + BsmHL/2], (14) 
dk' k 71 

where is the scattering angle of the neutrino in the laboratory system and is determined 
from the energies of the neutrinos in the initial(fc) and the final(/c') states as 

coseL = l- {k-k')Md/{kk'). (15) 

Since the deuteron is an iso-scalar object, only the iso-scalar vector current of the hadronic 
NC contributes to the ud elastic scattering. Neglecting a small contribution from the strange 
form factor, the matrix element of the hadronic current can be expressed with the iso-scalar 
elastic electromagnetic form factors of the deuteron. The structure functions, A and B, can 
be expressed with the electromagnetic Coulomb monopole (Gc), magnetic dipole (Gm) and 
quadrupole (Gq) form factors of the deuteron as 

B = ^Vil + v)G'm , (16) 
A = Gl + Iv'Gl + IvGl, , (17) 

with ?7 = Q'^/AMj and = (fc - k'f - (E^ - E'^y. We use the lA current to calculate the 
iso-scalar form factors. Since the dominant contribution to the elastic ud scattering is from 
Gc, effects of EXC is expected to be small in the energy region of our interest. 

In our numerical calculation, we use the ANLV18 potential j22| to generate the deuteron 
and two-nucleon scattering wave functions. The NN partial waves up to J = 6 are included 
for the deuteron breakup reactions. 



III. RESULTS AND DISCUSSIONS 

First we present the energy dependence of the total cross sections for the reactions in 
Eqs.(l)-(5). The total cross section for the neutrino and anti-neutrino deuteron reactions 
are shown in Figs. [1] (a) and (b), respectively. The neutrino CC reaction {yCC) gives the 
largest cross sections which are about 1/3 (1/2) of the neutrino-nucleon CC reaction at 
Ejj = 10(50)MeV. The pd elastic cross section is very small compared with that of uCC. 
In the high energy region around the pion production threshold {E^, ~ 140MeV), the pion 
production cross sections can be safely neglected j23|. 

An interesting quantity relevant to the supernova physics is a thermal average of the 
energy transfer cross section defined by[lH 

< (Tuj >T^ = J dEJ{n,E,)au{E,) . (18) 
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Ey [MeV] Ey [MeV] 

(a) (b) 

FIG. 1: Total cross sections for the neutrino-deuteron reactions. The sohd and dotted curves 
show the cross sections for (a) Ued — > e~pp {vCG) and vd upn (i^NC), respectively and (b) 
v'ed — > e'^nn (z^CC) and Pc? vpn (PNC), respectively. The dot-dashed curves show (a) the cross 
sections for v + n ^ e~ + p{vNCC) and (b) the cross sections for u + d ^ v + d {v — scatt). 



We assume here a Fermi-Dirac distribution for the neutrino, having in mind the neutrino 
flux from the supernova core. At the temperature with zero chemical potential, a neutrino 
with the energy E^, distributes as 

/(^^'^^) = (19) 

The energy transfer cross section auj{Ey) is evaluated by integrating the differential cross 
section multiplied by the energy loss of the incident neutrino with respect to the energy of 
the final lepton: 

auj{E,) = I dE[^E, , (20) 

for the CC reaction (absorption) and 

auj{E,) = J dK^{E, - El) , (21) 

for the NC reaction (scattering). This quantity is important for evaluating the neutrino 
heating behind the shock wave in the supernova. The energy of the incident neutrino 
is transferred to the matter (deuterons, nucleons, electrons/positrons) via the absorption 
(CC) or down-scattering {E^ < E^, NC). Note that the electrons, positrons and photons 
are regarded as a part of matter, being in thermal and chemical equilibrium with nucleons 
and nuclei in the supernova core, while the neutrinos are separately treated in the neutrino 
transfer calculations. 

The integrand of the average energy transfer cross sections, 

f{T,,E,)auj{E,) , (22) 
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FIG. 2: Energy [E^,) dependence of the energy transfer cross section for the CC vd reaction 
multiphed by the Fermi-Dirac distribution function at T;^ = 5 MeV(soUd curve) and 10 MeV(dashed 
curve) . 

for V CC is plotted in Fig. [2] as a function of the incident neutrino energy, E^, at = 5 
and lOMeV. The main contribution for the average cross section is from E,y ~ 20(60)MeV 
at = 5(10) MeV. We note that the contribution from the high energy tail of fau is 
appreciable. 

Finally our results on the thermal average of the energy transfer cross sections are shown 
in Fig. [3] and Table I. The largest cross section is due to i^CC which is almost the same 
magnitude as uCC The CC cross sections are an order of magnitude larger than the NC 
cross sections. We remark that the ud cross sections are appreciable even at relatively low 
temperature beyond the breakup threshold energy. This is a characteristic feature of the ud 
reaction involving the breakup of the weakly bound state with 2.2MeV binding energy. The 
rapid increase of the cross section at low energy makes the deuteron more preferable target 
than the other nuclear species in the supernova environment. 

In order to discuss the contribution of the deuteron in the neutrino heating mechanism 
through a comparison with other neutrino-nucleus processes, we evaluate the average energy 
transfer rate per nucleon through the neutrino and anti-neutrino reactions defined by 

< acj = ^[<(^^ >ucc + <(TLJ >pcc] , (23) 

< aUJ = ^[<(^^ >yNC + < CrtU >p7vc' + < O-u; >uscatt + <CrU >uscatt] , (24) 

where the factor 1/2 comes from the average over the neutrino and anti- neutrino. We show 
the calculated average energy transfer rate in Fig. |H For comparison we show the average 
energy transfer rate for A=3, 4 nuclei. For the CC reactions, the v-^Yi rate is calculated 
from i/-^H cross section given in Ref. [11 and the z/-^He rate is taken from Ref. [111 . For the 



NC reactions, we take the values for the z/-^He rate from Ref. ITj and the z/-^He rate from 
Ref. [l^ . It is remarkable that the average energy transfer cross sections for the vd reaction 
is significantly larger than those for ^H, ^He and ^He for the temperature range relevant 
to the supernova core. The vd rate is much larger at low temperature because of the low 
threshold energy of the deuteron as discussed previously. In the case of the CC reaction the 
vd rate dominates, being larger than those for A=3, 4 nuclei by orders of magnitude. 
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The large energy transfer rate of the deuteron at the wide range of the temperature sug- 
gests that the deuteron is potentially an important target for the neutrino heating mechanism 
in supernovae. The average energy transfer cross section for the CC reaction is sufficiently 
large to yield a significant heating rate for moderate neutrino temperatures. For example, 
the average energy transfer cross section for the CC reaction is larger than 10~^° cm^MeV 
at Tj, > 4 MeV, and gives a significant contribution to the heating rate in addition to the 
contribution of the neutrino-nucleon reactions. According to the expressions of the heating 



rate given by Haxton the heating rate due to the ud reactions for and Ue amounts to 
99.1 MeV/sec per nucleon, assuming T^^=Tp^=b MeV. The heating rate due to the vd reac- 
tions for /i and r type neutrino is 55.6 MeV/sec per nucleon at T^^=10 MeV. These values 
are 25-44% of the heating rate, 223 MeV/sec per nucleon, due to the z/-nucleon reactions 
at the same neutrino temperature of 5 MeV for v^./Vf.- In the above estimate, we assumed 
values of the neutrino luminosities (L,^=10^^ erg/sec) and the distance from the center (100 



km) in the expression for the heating rate given in Ref. [ll[. As the neutrino temperature 
is higher so does the heating rate, which makes the vd reactions more important. 

We remark that the deuteron as well as other light elements can appear in hot and dense 
matter in the heating region between the proto-neutron star and the shock wave. The Fe- 
group nuclei falling from the outer layer are dissociated into ^He and then nucleons at the 
shock wave. Through this dissociation, light elements [A < 3) including the deuteron can 



appear naturally in the hot environment |l7l . Il8| . Although the nucleons, having larger 



cross sections than those for nuclei, are the major targets of the neutrino after the complete 
dissociation, the deuteron and other light elements can contribute to the neutrino heating 
as an extra absorber of the neutrino. 

Light elements can be dominant targets in the case of enlarged shock radius due to 



the hydrodynamical SASI instability as shown by Ohnishi et al. [15| in multi-dimensional 
supernova simulations. In such case, the range of density and temperature is favorable to 
have abundant ^He with enough advection time to have neutrino absorptions. In their study, 
the i/-^He reactions may be crucial for the revival of shock wave in the marginal condition, 
but they require the enhanced u-^He reaction rates. Since the ud reaction rates is much 
larger than those of ^He as we have seen, the deuteron may play a similar role instead of 
^He during its dissociation. It would be interesting to explore the effect of ud reactions in 
the multi-dimensional supernova simulations with a detailed evaluation of mixture of light 
elements. The environment with an elongated shock wave due to the SASI instability may 
prefer to have abundant deuterons together with ^He. It has been shown that the deuteron 
appears (~ 1%) in the heating region behind the stalled shock wave in a realistic snapshot 
of central core after the bounce in the spherical supernova simulations 18 . 



Further studies on the suitable condition to have a sufficient amount of targets and a long 
enough reaction time periods are necessary to firm the possible role of neutrino-deuteron 
interactions for the success of supernova explosions. The systematic study of the abundance 
of light elements in hot and dense matter is now under way j23 |. 

In summary: the neutrino-deuteron reactions may have a potential impact on the neutrino 
heating mechanism in core-collapse supernovae. The cross section is large even at low 
neutrino energies relevant to supernovae because of the low breakup threshold for the bound 
state. The energy transfer rate through the i>d reactions is much larger than those of ^H, ^He 
and ^He even at low neutrino temperature. If deuterons appear abundantly in the heating 
region, as ^He in the SASI instability, they can contribute to the extra heating, which assists 
the revival of stalled shock wave. It is interesting to explore whether the ud heating can 
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Tv [MeV] Tv [MeV] 

(a) (b) 

FIG. 3: Thermal average of energy transfer cross sections. The sohd and dotted curves in (a) 

((b)) show the cross sections for Ued e~pp{iyCC) and i^d z/p?i(i/NC) (Pgd e+nn(z/CC) and 

i/d — > Pj>n(^'NC)), respectively. The dot-dashed curve in (a) and (b) shows cross section for the 
elastic I'd scattering. 




FIG. 4: Averaged energy transfer cross sections in unit of 10~^^ MeV cm^. The solid, dashed and 
dash-dotted curves show the ud, z^-^He and i^-'^H (left panel) and z/-^He (right panel) cross sections, 
respectively. (See the main text for the references on A=3, 4 nuclei cross sections.) 



affect the dynamics of shock wave in reahstic supernova simulations which take into account 
the mixture of light elements and the associated neutrino reactions. 

The evaluated data of the ud energy transfer rate will be available on the web site for 
the neutrino deuteron reactions 25|]. The data can be used to implement these processes 
into the heating rate in supernova simulations. Detailed neutrino differential cross sections 
based on the current theoretical approach can be provided for numerical simulations of the 
neutrino transfer which require the angle and energy variations. It will be also interesting 
to study the lyd processes around the surface of proto-neutron stars and its influence on the 
neutrino nucleosynthesis. The other channels of weak processes on the deuteron are now 
under investigation. 
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T^(MeV) 



V CC 



V NC 



V NC 



V Scatt. 



5.11 X 10"^ 
7.95 X 10-3 
3.89 X 10-2 
1.18 X 10"^ 
2.75 X 10"^ 
5.42 X 10"^ 
9.51 X 10-1 
1.53 X 10° 
2.31 X 10° 
3.30 X 10° 



1 
2 
3 
4 
5 
6 
7 
8 
9 

10 



1.06 X 10° 
1.51 X 10^ 
6.37 X 10^ 
1.71 X 10^ 
3.67 X 10^ 
6.79 X 10^ 
1.14 X 10^ 
1.78 X 10^ 
2.64 X 10^ 
3.73 X 10^ 



2.68 X IQ-i 
7.50 X 10° 
3.63 X 10^ 
1.01 X 10^ 
2.14 X 10^ 
3.87 X 10^ 



1.43 X 10-1 
1.66 X 10° 
5.68 X 10° 
1.31 X IQi 
2.50 X IQi 
4.28 X IQi 



6.29 X 10^ 6.80 X 10 



9.49 X 10^ 
1.35 X 10^ 
1.85 X 10^ 



1.02 X 10^ 
1.48 X 10^ 
2.06 X 10^ 



1.38 X 10-1 
1.56 X 10° 
5.17 X 10° 
1.15 X lOi 
2.13 X lOi 
3.51 X lOi 
5.37 X lOi 
7.80 X lOi 
1.09 X 10^ 
1.46 X 10^ 



TABLE I: Averaged energy transfer cross sections < ctoj > in unit of 10 MeV cm^. 
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